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Abstract 

Among various super symmetric (SUSY) standard models, the gravity mediation 
model with a neutralino LSP and the gauge mediation model with a very light grav- 
itino are attractive from the cosmological view point. These models have different 
scales of SUSY breaking and their underlying physics in high energy is quite dif- 
ferent. However, if the sparticles' decay into the gravitino is prompt in the latter 
case, their collider signatures can be similar: multiple jets and missing transverse 
momentum. In this paper, we study the discrimination between these models in 
minimal cases at the LHC based on the method using the significance variables in 
several different modes and show the discrimination is possible at a very early stage 
after the discovery. 



1 Introduction 



Supersymmetric (SUSY) standard model (SM) is the most promising candidate for the 
physics beyond the standard model (BSM). Once its presence is inferred from experiments, 
the next problem is to search for the underlying mechanism of SUSY breaking. The scale 
of SUSY breaking is one of the most important clue for this search. 

There are many proposed models with various SUSY breaking scales. A major candi- 
date is the gravity mediation model with a medium SUSY breaking scale. In the model, 
the lightest SUSY particle (LSP) is the Bino-like neutralino. It is stable and becomes a 
good candidate of a WIMP dark matter (DM), if i?-parity conservation assumed. At the 
LHC, the LSP is observed as a missing transverse momentum (pt)- The missing serves 
as a distinctive signature as well as high-pr jets and other decay products coming from 
decay chains of produced SUSY particles. 

Another attractive candidate is the gauge mediation (GM) model with a low SUSY 
breaking scale since it solves the SUSY flavor problem. The gravitino becomes the LSP 
in the GM model and the collider signatures strongly depend on the gravitino mass 
m3/2, which is related to the SUSY breaking scale a/F as ^3/2 = F / {-s/^Mpi), and the 
nature of the next LSP (NLSP). In the case of the heavier gravitino, the SUSY event 
at the LHC involves a non-pointing NLSP decay [Tl, ^2] or a long-lived NLSP [3J, which 
becomes a characteristic signature. On the other hand, if the gravitino is very light, 
the NLSP decay is prompt and the gravitino is observed as a missing pt- Then, the 
event signature resembles that of the gravity mediation model. From the cosmological 
point of view, such a very light gravitino with 7713/2 ^ 10 eV is very attractive since 
the model is free from gravitino problems So the discrimination between the gravity 
mediation model and the GM model with a very light gravitino is an important problem. 
However, this discrimination is difficult. In principle, measuring the LSP mass provides 
good information for discrimination. But, to measure the LSP mass precisely, a large 
integrated luminosity is required, and in some case, obtained result may have degenerated 
solutions. 

In this paper, we concentrate on the minimal types of the gravity mediation model 
and the GM model. We call these models the mSUGRA model and mGMSB model. 
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respectively. The NLSP in the inGMSB model becomes either the Bino-like neutralino 
Xi or the lighter slepton i (stau fi). In the former case photons from the Bino decays are 
produced and in the latter case leptons are produced, and the photons or leptons are ex- 
pected to be the discriminant of the GM model. High energy photons are rarely produced 
in the gravity mediation model and there is no difficulty in discriminating the GM model 
with a Bino NLSP. However, leptons are also produced in the gravity mediation models 
and thus the discrimination has a more quantitative nature. Also, in the experiment, the 
signal events are smeared by the SM backgrounds and the effect of detector efficiencies. 
Furthermore, if the lightest neutralino mass rests between the lighter slepton and stau: 
rriri < ^ ^Ir' lepton emissions can be suppressed and/or hidden. 

The purpose of this paper is to study this discrimination between the mSUGRA and 
the mGMSB models in the realistic experimental setup. We propose a simple method 
of discrimination and show the integrated-luminosity dependence of the degeneracy of 
models. The method is applicable at an early stage of the LHC soon after the discovery. 

The organization of the paper is as follows: In the next section, we describe the 
parameter spaces and the mass spectrums of the models considered in this paper and 
introduce the collider signatures of these models. We present the result of our Monte 
Carlo simulation and show the idea of using significance variables of multiple modes in 
Section El In Section HI we discuss the discrimination of the mGMSB and the mSUGRA 
models based on a test on significance variables. Section |5] is devoted for discussions and 
conclusions. 

2 mSUGRA vs mGMSB: event signatures 

In this section, we first review the parameter space of each model and discuss the LHC 
signatures of these models. Then we describe the idea of discriminating these models at 
an early stage at the LHC. 

2.1 Model parameters and mass spectrum 

In this paper, we compare the collider signatures of the mSUGRA model and the mGMSB 
model. First, we describe the model parameters of each model. 
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The mSUGRA model represents a mass spectrum of the gravity-mediated SUSY break- 
ing model assuming the flavor blindness for the scalar masses and the CP property. The 
model parameters are mg, "011/2, tan /3 and sign /i, which is the unified scalar mass, the 
unified gaugino mass, the universal A-term, the ratio of the Higgs VEVs and the sign of fi, 
respectively. The low-energy mass spectrum is obtained by solving the renormalizat ion- 
group equations with these input values at the GUT scale. The resulting LSP is typically 
the lightest neutralino or the lightest slepton. From the cosmological view point, we only 
consider the case of a neutralino LSP. 

As a parameter space of the mCMSB model, we adopt a parametrization used in 
Ref. [5]. Here, the ratio between the gaugino mass and the sfermion mass is left free while 
the "GUT relations" among the gaugino masses and the sfermion masses are imposed. The 
A-term is set to be and the S-term is tuned to give a correct low-energy electroweak 
parameters. The model parameters are A^, Ag, Mmes, tan/3 and sign /i. The sfermion 
masses and the gaugino masses at the messenger mass scale Mmes are given by 

M. = (1) 

a 

where (a = 1,2,3) are the SM gauge couplings and Ca{i) are the quadratic Casimir 
coefficients for each field. The A^ and A^ are related to the parameters of the conventional 
mGMSB model parameters as 

A^ = 2A^5A^ (3) 
A, = N,A, (4) 

where is the number of messengers in the 5 + 5* representation and A = F/M^es is 
the ratio of the square of the SUSY-breaking scale a/F and Mmes, when F/M^^^ ^ 1. 
So in other words, the parameter representation in Eqs. ([1]) and (|2]) is the same as the 
conventional one with an arbitrary (non-integer) number of messengers. 

2.2 mSUGRA signal 

The LSP is the lightest neutralino in the mSUGRA model as described in the previous 
section. Once the SUSY particle, mostly the gluino and the squarks, is produced it 
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decays in a cascade down to the SM particles and the LSP. Since the stabihty of the LSP 
is guaranteed by the -R-parity conservation, the produced LSPs then escape outside the 
detector and are observed as a missing p^- 

Since the gluinos and the squarks are colored particles and the LSP is a color-neutral 
particle, the cascade decay products include jets, as well as other SM particles depending 
on the detailed form of the decay chain. There are many possible decay chains. For 
example, in the decay chain — )■ Xol ~^ Xo^^Qy leptons are emitted as well as jets and 
missing p^. Likewise, the decay products may include leptons, tau-jets or 6-jets as well 
as the generic signal of jets and missing momentum. 



In the mGMSB model we are interested in, the LSP is the gravitino. As in the mSUGRA 
case, the produced SUSY particles decay into the gravitino LSP. In the mGMSB case, 
however, the role of the NLSP becomes more important. 

The decay length of the NLSP into the gravitino LSP is given by 



Since this is much longer than the one resulting from the SM interaction, a produced 
SUSY particle usually decays into the NLSP before it directly decays into the LSP. 
Thus, the event signature of the mGMSB model strongly depends on which particle is 
the NLSP. 

As mentioned in the previous section, the NLSP in the mGMSB model is either the 
lightest neutralino or the lightest slepton (stau). In the neutralino-NLSP case the NLSP 
emits a photon and in the slepton-NLSP case the NLSP emits a lepton with a sizable 
branching fraction. The typical mGMSB signal, therefore, contains multi-jets + missing 
-|- multi-photons or multi-jets + missing + multi-leptons. 

The lightest slepton is the lighter stau fi and the NLSP in the slepton-NLSP case. So 
as well as the multi-leptons signal, we expect events with multi tau-jets. The multi-tau-jets 

^ There are exceptional cases. One example is the case when the right-handed slepton mass is degen- 
erated to the lighter stau mass for small tan /3 [S] . Another example is the case that the masses of the 
lightest neutralino and the lighter stau are degenerated for Ag ^ A^. 



2.3 mGMSB signal 




(5) 
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mode becomes more important in the large tan /3 case, since then the heavier sparticles' 
decays into staus becomes more dominant because of the kinematical preference and the 
enlarged coupling of the fi due to the large left-right mixing of the staus. 

2.4 Strategy for discrimination 

As described above, an important difference between the mSUGRA model and the mGMSB 
models arises from the important role of the NLSP decay in the mGMSB model. In 
the mGMSB model with a neutralino-NLSP, the signal event usually accompanies two 
photons, and in the case with a slepton-NLSP, the signal accompanies multi-leptons. 
Grounded on these facts, we compare the models by using the following modes: no lepton 
+ 2 jets, no lepton + 4 jets, same-sign (SS) 2 leptons, 2 tau-jets and 2 photons modes. 
In table [H we summarize the naively expected number of signal events in these modes 

Table 1: Naively expected number of signal events. The mark @ indicates that a large 
number of signals are expected and the mode is qualitatively "good" for model search. 
Similarly, Q, A and x indicates that the corresponding mode is qualitatively "not bad", 
"not good" and "bad" . 



mode 


mSUGRA 


mGMSB f-NLSP 


mGMSB x?-NLSP 


Qi + (2 jets/4 jets) 


® 


A 


A 


SS 2i 





® 


A 


2t 


A 


® 


A 


27 


X 


X 


® 



for each model. For example, in the mSUGRA model, the no lepton + (2 jets/4 jets) 
mode is expected to be the most viable and the modes of two leptons or tau-jets are 
the next, but the two photons mode is bad because photons are hardly expected in this 
model. Similarly, the two leptons mode or the two tau-jets mode is expected to be the 
most viable one for the mGMSB model with a slepton-NLSP and the two photons mode 
is the best mode for the mGMSB model with a neutralino-NLSP. Since the location of 
the viable mode differs for the three cases, we may distinguish these models by observing 
which of the four modes shows the most signal events or which does not. 

However, we need sophistication to apply this idea in the realistic experimental set- 



6 



up. First, we must include contamination by the SM backgrounds and the detection 
efficiencies. The signal events may be hidden unless we carefully reduce the SM back- 
ground contribution, which might dilute the apparent differences of characteristic events. 
Further, there are a certain amount of misidentification and fake rate of particles by the 
detection systems. This effect is non-negligible in the case of many leptons and photons, 
and is also a source of dilution of the signal events. Second, we need a quantitative way to 
describe the distinctions among various models, which can incorporate the experimental 
errors. 

In the next section, we perform a Monte Carlo simulation, which incorporates the SM 
backgrounds and the detection efficiencies. We use the significance variable (see Sec. 13.21 
for the definition and explanation) to solve the above problems. 

3 mSUGRA vs mGMSB: simulation 

In this section, we ffist describe the set-up of our Monte Carlo simulation and introduce the 
significance variable. We also review the current experimental constraints. We then show 
the result of our simulation and discuss the possibility of discrimination with significances 
in several different modes. 

3.1 SM backgrounds and detection efficiencies 

For the multi-jets mode, the main SM backgrounds come from the QCD jets, tt and 
gauge boson production events. For the multi-lepton mode tt and W+ jets are the main 
background and for the multi-photon mode, the 77 events are. To estimate these back- 
grounds and other sub-leading backgrounds, we used the programs MC@NLO 3.42 [7] (for 
tt, WW, WZ and ZZ), Alpgen 2.130 (for Wj, Zj and W/Z + bb/tt), MadGraph 4.1.44 
[9] (for ti/W/Z + -f/-f-f), and Pythia 6.4 ^ (for 77 and QCD jets). As for the detec- 
tor simulation, we used the AcerDet 1.0 [11]. For the isolation criteria for leptons and 
photons, we used the default setting of AcerDet. 

In our simulation, we are interested in the signals with leptons or photons. The 
detection efficiency and misidentification of these particles are important in estimating 
both the signal events and background events. In the following simulation, we include the 
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misidentification and fake rates of the leptons, photons and jets. The rates we have used 
in the analysis are shown in Table |2l The values are taken from Ref . |20] . For the SM 
backgrounds, however, the efficiencies are chosen to give more conservative results than 
the actual experimental data. 



Table 2: Misidentification and fake rates. The efficiencies are taken from Ref. [2^. For 
SM backgrounds, the efficiencies are chosen to give a conservative estimation. 





j — 7- e 


J ^ 


J ^ 7 


T — ^ j 


j 


e — 7- j 
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SUSY 


0.3 % 


0.3 % 


0.02 % 


60 % 


1 % 


27% 


30 % 


3 % 


20 % 


SM BKG 


0.3 % 


0.3 % 


0.02 % 


20 % 


1 % 


% 


% 


3 % 


% 



3.2 Significance 

Under a fixed set of event cuts, we estimate the event number of the signals Ng and the 
backgrounds A^;,. To compare many model points in several modes on the same footing, we 
use the significance variable Z, which is a function of A^^ and A^;,. We use the significance 
incorporating both the statistical and systematic uncertainties of the backgrounds and 
adopt the one used in Ref. [20]. The systematic uncertainties of the backgrounds are 
taken to be 50% for the QCD multi-jet background and 20% for others. 

Given the number of the signal events A^^ and the background events A";, with the 
uncertainty 5Nb, the significance is given by calculating the convolution of the Poisson 
distribution with some "posterior" distribution function. As the posterior distribution, 
we take the gamma distribution as suggested in Ref. [12j. The resulting significance Z is 
given by |12] 

Z = V2erf"^(l -2j9), (6) 

with 

_ B{N, + AT,, 1 + NI/5NI 5Nl/{N, + dN^)) 
^ B{N, + AT,, 1 + N^J5N^,) ' ^ > 

where erf~^ is the inverse error function and 

B{a,b,x)= [ dtt^-\l-tf-^ (8) 
Jo 
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is the incomplete beta function and B{a, b) = B{a, b, 1) is the usual beta function. If we 
take the limit 6Nb — )■ 0, the Eq. ([7]) reduces to the probability in the Poisson distribution. 
The statistical error of the significance SZb is depends on both Ng and Nf,. The error is 
(9(1) for the moderate values which we are interested in. 

3.3 Experimental constraints 

Up to now, the most stringent experimental constraints are given by the ATLAS and CMS 
experiments at integrated luminosities about 1 fb^^. Analyses relevant for the GMSB and 
mSUGRA models are searches of jets with no lepton [13] or one lepton [14] by the ATLAS 
collaboration and searches of same-sign [15] or opposite-sign ^Jj dilepton and photons 
[T7] by the CMS collaboration. 

We show the constraints given by these experimental data in Fig. [1] To obtain the 
constraints, we apply the 95% C.L. exclusion limits given by the above analyses to our 
model points. To incorporate the NLO corrections to the cross section, we multiply a 
uniform /c-factor of 1.5 for all model points, except for the mCMSB model points with a 
Bino-NLSP, at which the dominant contribution to the SUSY production is given by the 
direct Wino production. 

3.4 Signal events 

For both the mSUGRA and mCMSB models, the MSSM mass spectrums and the decay 
tables are calculated by using the program ISAJET 7.72 [T8]. The SUSY events are 
generated by using the Herwig 6.510 [19] and we have used the AcerDet 1.0 [Hj for the 
detector simulation. 

For the mSUGRA model, the grid in parameter space is taken to be < mo < 1200 
GeV (in steps of 20 GeV), 100 < mi/2 < 600 GeV (in steps of 20 GeV), = 0, mo, -mo, 
tan/3 = 10, 40 and sign /i = +. For the mGMSB model, the grid in parameter space is 
-(100)2 < A2 < (400)2 TeV2 (in steps of 5 TeV), (100)^ < A^ < (400)^ TeV^ (in steps of 
5 TeV) and tan (3 = 10, 40. The messenger mass scale is fixed as mmes = 200 TeV 
the gravitino mass is set 1713/2 = 16 eV. 



and 



^ To avoid a messenger tachyon, the condition F < mf^^^^ is required, which indicates mmcss ^ 100 
TeV. 
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Figure 1: Experimental constraints given by the LHC experiment at an integrated lumi- 
nosity of about 1 fb~^ at the 7 TeV run. We also show the benchmark points PTl-10. 
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Because we are interested in the model points that are potentially discovered at an 
early stage of the LHC and are not excluded by the experiments, we select model points 
within a certain mass region to reduce the size of computational calculation. For the 
mSUGRA, we select model points with 800 < mj^ < 1600 GeV and 400 < m-g < 1400 
GeV. For the mGMSB, we select model points with 850 < m^^ < 1600 GeV for the 
slepton-NLSP case and 220 < m^o < 350 GeV for the neutralino-NLSP case. The number 
of model points after this selection is 5633 (4046) for the mSUGRA (mGMSB) model. 
After imposing the experimental constraints reviewed in the previous section the number 
reduces to 3654 (3048). 

3.5 Modes and kinematical cuts 

In the following analysis, we use five different modes, no lepton modes with two and four 
jets, two photons mode, two r-jets mode, and the same-sign (SS) two leptons mode. 
The kinematical cuts for these search modes are summarized in Table |3l In the table. 

Table 3: List of kinematical cuts used in the analysis. The i?™'^^ is the absolute value of 
the missing px and Meg is defined in Eq. (jO]). We require isolation of leptons and photons 
(see text). In the four jets with no lepton mode, the values in brackets are chosen for each 
model point to optimize the significance. 



mode 




SS 2£ 


At least two leptons with > 20 GeV. The leading two leptons have 
the same charge. At least 2 jets with px > 50 GeV and for the leading 
jet px > 100 GeV. Ef''"" > 100 GeV, Mes > 800 GeV. 


2r 


At least two tau-jets with px > 40 GeV. At least 4 jets with 
px > 50 GeV and for the leading jet px > 100 GeV. 
E^'^^ > 100 GeV, Meff > 1000 GeV. 


27 


At least two photons with px > 60 GeV and for the leading photon 
Px > 90 GeV. EJ^'^^ > 100 GeV 


0£ + 4 jets 


At least four jets, requiring px > {50, 100} GeV, and for the leading jet 
Px > 100 GeV. No leptons with px > 10 GeV. E™'^^ > 200 GeV, 
Meff > {1000, 1200} GeV and Ef^^^M^s > {0.2,0.3}. 


0£ + 2 jets 


At least two jets with px > 200 GeV, and for the leading jet 
Px > 300 GeV. No leptons with px > 10 GeV. E^'^^ > 400 GeV, 
Meff > 1200 GeV and Ef^^^Mes > 0.35. 
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the E^^^^ denotes the absolute value of the missing px, and the effective mass is Meg 
defined by 

2 or 4 

Meff = i^r^ + ^ ptU) + J2 pt(^), (9) 

jr: jet £: lepton 

where the summation on the jet pt is for the leading two jets in the two r-jets mode and 
the no lepton with two jets mode, and for the leading four jets in other modes. The cuts 
in each mode are those typically used by the current experiment at the LHC. To fix the 
numerical values for kinematical variables, we performed simple optimization to maximize 
the significance. In the four jets with no lepton mode, for which the maximized values 
depend rather heavily on the location of the model points, we left several cuts unfixed. 
When a cut has two values in brackets, it means that the kinematical cut is optimized 
within the described values. 



3.6 Example model points 

We first pick up several example model points to illustrate the result. The 10 selected 
model points (PT 1 to PT 10) are depicted in Fig. [1] and the corresponding model pa- 
rameters are summarized in Table |H 



Table 4: Model parameters of the example model points. 



mGMSB 


As (TeV) 


Ag (TeV) 


tan /3 




PT 1 


350 


200 


10 




PT 2 


350 


200 


40 




PT 3 


150 


180 


10 




PT 4 


150 


120 


40 




PT 5 


100 


260 


10 




PT 6 


120 


230 


40 




mSUGRA 


mo (GeV) 


mi/2 (GeV) 


tan/3 




PT 7 


260 


500 


10 





PT 8 


260 


500 


40 





PT 9 


1100 


260 


10 





PT 10 


1100 


260 


40 






PT 1 to 6 are mGMSB model points. PT 1 and 2 represent models with a Bino-NLSP, 
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Figure 2: Significances Z's of tlie example model points for the five modes listed in Tab. |3l 
The size of red boxes shows the corresponding significance with an integrated luminosity 
of 30 fb"^ at 7 TeV run. 



PT 3 and 4 models with a slepton-NLSP but the fi mass and the lightest neutralino mass 
are nearly degenerate, and PT 5 and 6 with a slepton-NLSP without such a degeneracy. 
PT 7 to 10 are mSUGRA model points. Here, we take Aq = 0. PT 7 and 8 are typical 
models with small mo and large mi/2, and PT 9 and 10 are typical models with large mo 
and small mi/2- In the Appendix, we list the mass spectrum of relevant superparticles 
and their dominant decay modes in Tab. [TOjandfTTl 

In Fig. [21 we illustrate the significances of the example model points in the five mode 
in Tab. [3] with the integrated luminosity of 30 fb~^ at 7 TeV run. The number of signal 
events and our estimation of the SM background is listed in Tab. [T2]in the Appendix. We 
see that the naive expectation in Table [H is in fact visible in terms of the significances. 
An exception is PT 4, where the no lepton modes have unexpectedly large significances. 
At this point, the lighter stau mass is nearly degenerated to the lightest neutralino mass 
and the significant proportion of taus or leptons from the upper-stream of decay-chains 
become soft. This fact results in smaller events with high-px leptons. 
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3.7 Results for all model points 

We now illustrate results for all model points. Since the significance variable is an increas- 
ing function of the number of signal events after the cut, its absolute value scales with the 
total cross section of the model point. Thus the meaningful quantity for discrimination 
is the ratio of significances in different modes or a similar quantity. In the following we 
show the distribution of the significance in each mode in comparison with the one in the 
four jets with no lepton mode Z(0£4j). The significances are calculated for the integrated 
luminosity of 30 fb^^ at 7 TeV run. 

We first compare the photon mode significance. Since the multi-photon signal is a 
very distinctive signature for the mGMSB model with a neutralino-NLSP, the model 
is expected to be distinguished clearly from the others in terms of the photon mode 



significance. In Fig. 3(a) , we plot the photon mode significance Z{2'y) versus the four jets 
mode significance Z{0i4j). Each point in the figure corresponds to a model point. A red 
and triangle (blue and circle, black and square) mark represents a mSUGRA (mGMSB 
with a slepton-NLSP, mGMSB with a neutralino-NLSP) model point. We can check a 
clear separation of the mGMSB model points with a neutralino-NLSP from the other 
models as expected. A few model points of the GMSB model with a slepton-NLSP have 
the significance larger than five. At these model points, the lightest neutralino Xi and the 
lighter slepton are almost degenerated and the decay of the lightest neutralino Xi into the 
lighter slepton is kinematically suppressed. Thus, these model points are parametrically 
very near to the neutralino-NLSP case and we here treat them as the neutralino-NLSP 
case. 

Next we show the significance in the SS two leptons mode Z{SS2i). To clarify the illus- 
tration, we deselect the mGMSB model points with a neutralino-NLSP (or more precisely 



those with Z(2'y) > 3). In Fig. 3(b) (and the following figures), the rest mGMSB model 
points are marked with blue and circle and the mSUGRA model points are marked with 
red and triangle. We see a clear inclination that the mSUGRA model points have smaller 
Z(SS2£)/Z(0£4j) and the mGMSB model points have larger Z {SS2i) / Z (OiAj) , showing 
that for most of the model points the lepton richness is indeed a good discriminator even 
with the smearing due to the SM background and detector effects. However, we also find 
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Figure 3: Scatter plots of significances Z's for an integrated luminosity of 30 fb^ at 7 
TeV run. The black square points are the mGMSB model points with a neutralino-NLSP, 
the blue circle points the mGMSB model points with a slepton-NLSP, the red triangle 
points the mSUGRA model points. 
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many mGMSB model points with rather small Z {SS2i) / Z {0i4j) . As discussed in the 
previous section, a typical reason for a smaller Z{SS2i) is nearly degenerated NLSPs. If 
the slepton-NLSP is nearly degenerated to the lightest neutralino, then a possible leptons 
coming from the decay of the lightest neutralino into the slepton becomes soft and difficult 
to be detected. Another but related cause for a less lepton signature is larger tan/3. As 
indicated in the previous section, when the tan /3 is large, the SUSY decay chains involve 
more staus than other sleptons, because large tan /3 induces the large left-right mixing of 
the staus, which leads to the small stau mass and large interaction originated from the 
component of the left-handed stau. Thus decay into the stau is preferred from both rea- 
sons of kinematics and interaction strength. Note that the decrease of the leptons results 
in not only a smaller Z{SS2i) but also a larger Z{0i4j), showing that a smaller Z(SS2i) 
at those points is not cased by a mere scale of cross section. 

To rescue the mGMSB model points with a small lepton mode significance, we now 



show the significance in the two r-jets mode versus that in the four jets mode in Fig. 3(c) 



Here we see an apparently better separation between the mSUGRA and the mGMSB 
models compared to the dilepton mode. Especially, those mGMSB model points with 
larger Z{0i4j) have larger Z{2t). This is another check that in those model points the 
leptons are replaced by r-jets. For the discrimination, this r-jets mode is important 
since it plays a complementary role to the dilepton mode. Note that for most of the 
mGMSB model points, the Z{2t) is small compared to Z(SS2£). This is because of the 
low detection efficiency of the r-jet. 



Finally we look at the two jets mode in Fig. 3(d) The mSUGRA model points are 
spread over a wide region, but the mGMSB model points have smaller Z{0i2j)/Z{0i4:j). 
The two jets mode alone does not work as a discriminator, but it strengthen the separation 
between the mGMSB model points and part of the mSUGRA model points. In summary, 
with the inclusion of the SM background and detection efficiencies, the photon mode, SS 
dilepton mode together with two r-jets mode and their counters, two jets and four jets 
with no lepton modes, remain serving as discriminator for the mSUGRA and mGMSB 



models. In Fig. |3(e)]we plot two sums Z(SS2£) + Z{2t) vs. Z{OUj) + Z{0i2j) for both 
models. We see a fairly good separation between the two models. 

In reality, the significance values have uncertainties due to the experimental errors and 
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errors from the MC statistics. With an increasing integrated luminosity, each model point 
gains larger absolute values of the significances and thus the distance from the opposite set 
of models (mSUGRA vs. mGMSB) become larger. Then the next problem is how much 
the integrated luminosity is needed in order that the model point is well-separated from 
the opposite models. In the next section, we study the integrated-luminosity dependence 
of the separation between the two model. 

4 Significance test for discriminating models 

Given the above results, we now go on to the discrimination of the mGMSB and the 
mSUGRA models. Discrimination here means that excluding the mGMSB models from 
the mSUGRA models and vice versa. In the previous section, we have seen that even with 
the inclusion of the effects of the SM background events and the statistical uncertainty, we 
expect a fairly good separation between the two models at earlier stages of the experiment 
after the discovery. However, to answer the question of how far and at which luminosity 
the discrimination is possible we need a further machinery. In this section, we develop a 
statistical method and answer this question. 

The statistical method we use is based on the test on the significance variables. 
Given two model points, rria and m^, and an integrated luminosity L, the variable is 
defined as 



where i runs for modes used in the analysis and 6Z is the error of Z. In general, 6Z 
contains various systematic errors as well as the statistical error. In the following analysis, 
we only consider the expected experimental statistical error and the uncertainty from MC 
statistics. We may also define by using functions of significance variables. Then the 
denominator of Eq. ( ITOj) is replaced by the error of corresponding functions. 

If, for the two model points, the calculated is smaller than a certain value, their 
difference of significances are considered as statistically insignificant. In the following 
analysis, we use a p- value of 95% for this criterion. The corresponding value depends on 
the number of degrees of freedom, i.e. the number of significances used for the calculation 




(10) 
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of x^. We say that two model points are in the vicinity of each other if their is smaller 
than than this value. 

Now, we consider the problem of testing whether a certain model is consistent 
as an mGMSB model or an mSUGRA model. We first prepare a set of model points in 
these models. We denote by M the set of model points and indicate as M = G the set of 
mGMSB model points and M = S that of the mSUGRA models, respectively. Then let 
us perform a test and collect those model points of M in the vicinity of ttiq. We denote 
this set of model points as Suifno). If Suifno) is not empty, there is a certain probability 
that mo is in the model of M. However, we must take account of the density of points 
in M. If there are a large number of points in M, there might be a high probability of 
the 5'm(?«o) being non-empty even if mo is not a good candidate in the model of M. On 
the other hand, if the cardinality of M is small there may be a possibility that 5'm(^«o) 
is empty even if mo is actually in the model of M . 

For this purpose, we first consider a point m in M as mo. For the care of this last 
problem, we must prepare a dense enough set M . Then for each point m in M, SMiiT^) 
contains at least a few points. The number of points in SMijn) indicates the sensitivity 
of the significances in the vicinity of m. We now consider a quantity sm{i^) defined as 



If the variation of ^j^StyiiiT^) is smooth compared to the density of model points around 
m in the region of Suijn)-! the value of SM{rn) is near unity. Similarly, we define for a 
general model point mo 



If the density of M is high enough and the variation of #<S'm(?77.) is smooth, we expect 
that sm(?71o) is also near unity if mo is indeed in the model of M. Note that this con- 
dition is stronger than the mere condition that SMiiTio) is not empty. In this method of 
discrimination, the condition for mo to be consistent with being in the model of M is that 
sm("^o) is near to unity within the errors, which are discussed further below. 

Before discussing the errors, we have several remarks on this method. Let us consider 
the limit of infinite density of M and infinite luminosity. Then the distribution of model 




(11) 




(12) 
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points in M in the significance space spreads over a certain region and then sm for a 
given model is unity or zero, indicating the model is in or out of this region. In other 
words, this method is a sort of in-or-out test and does not take the density distribution 
of M in this limit. So this method is useful to discriminate well separated models, but 
may be useless for those models, in which a tiny amount of model points spreads over 
a wide region. In practice, we have a finite density of M and a finite luminosity, and 
the region is approximated by discrete points. In this case, the sufficiency of the density 
of M may be tested by SM{fn) for m in M as noted above. There may be some points 
with a small S'm(^) compared to the rest. One reason for such a case is the existence of 
points, called "cliffs" in Ref. [21]. Around these points the variation of the significance 
is large compared to that of the parameter. The boundaries in the parameter space or 
in the mass spectrum space may be a cause for such a case. This case is often physically 
interesting and may require further and detailed studies. 

The errors come from the effects of both the finite density and luminosity. In a case 
of finite luminosity, we have a certain variation of the value of the sm('^o) other than 
unity or zero if ttiq is near the boundary of M in the significance space. Assuming M is 
uniformly distributed, the sm(^o) at the boundary is equal to log 2 ~ 0.69 (independent 
of the dimension of the significances) and increases as rriQ moves from the boundary into 
the region of M to a maximal of ~ 1.2 depending on the dimension. Away from the 
boundary sm(^o) decreases to zero. These values may vary well according to the shape 
and distribution of concerning model M. Since we do not know the distribution of M 
and the position of mo in the significance space a priori, we handle this variation as an 
uncertainty and we take the value (5;Sm(^o) = 0.5 for this uncertainty in the following 
analysis. There are also errors on sm('^o) coming from the finite density of M. We 
adopt a rough approximation of this error by assuming that S'A/(mo) shapes like a multi- 
dimensional sphere in the space of significances. In this approximation, the error of 
SA/(mo) is given by 




(13) 



where d is the degrees of freedom of and Cd = vr' 



'^/Vrd + 1)- As a result, the condition 
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for the consistency of mo in the model M can be written as 

kM("^o) - 1| < SsMimo) = 6gSM{moy + 5/SA/(mo)2. (14) 

We remark that this estimation of 6s m is very rough and we assume that the variation of 
Sm in the parameter space is smooth. 

We examine this method of model discrimination with the result of the Monte Carlo 
simulation in Sec. 13.71 We use as the set M = G or 5* the set of all simulated model points 
described in Sec. 13.41 As test model points, we choose from the same set model points 
which pass the experimental constraints and satisfy the discovery condition, for which we 
require at least one of the significances is greater than five for 30 fb^^. The set of test 
model points are denoted by G' and S'. The number of model points are 4046, 1091, 5633 
and 754 for G, G', S and 5", respectively. We show the results for calculated in three 
different combinations of the modes listed in Table El The first one uses the lepton 4 
jets mode, 2 photons mode and SS 2 leptons mode, the second one uses the sum of two 
lepton modes, 2 photons mode and the sum of 2 r-jets and 2 leptons modes, and the 
third one uses all the five modes. We call these three cases as Tl to T3. The results are 
illustrated with integrated luminosities L = 1, 3, 10, 30 and 100 fb~^ in the 7 TeV run. 

We first check model points with small #S'm- In Tabled the number of model points 
in M = G and S with ^Sm less than ten is presented. Because most of the points with 

Table 5: The number of model points with < 10, which pass the experimental 

constraints. The total number of model points are 3048 for G and 3654 for S. The 
numbers for L = 1 and 3 fb~^ are all zero. 



L 


10 fb-i 


30 fb-i 


100 fb-i 


G in Tl 





1 


2 


G in T2 


1 


3 


3 


G in T3 





2 


4 


S in Tl 











S in T2 











S in T3 












"^Sm < 10 are near the artificial boundary in the model parameter space, especially in the 
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low mass region, and they are unlikely to affect the result of the test, only those which pass 
the experimental constraints are counted in the table. The few mGMSB model points in 
the table are those nearly on the boundary where the NLSP changes. The vicinity of these 
points is a good example of a "cliff". The followed result shows the effect of this region 
for discriminating the models, but a further study might be important if the experiment 
reveals its possibility. 

Now we present our result of the tests Tl to T3. In Table [6], sg and ss averaged in 
G' or S' are presented. First, we see that sm calculated in the corresponding model 

Table 6: Evolution of the averaged sg,s with respect to the integrated luminosity L for 
significance test Tl to T3. 





L 


1 fb-i 


3 fb-i 


10 fb^i 


30 fb-i 


100 fb-i 


Tl 


Sg averaged in G' 


1.036 


1.042 


1.011 


1.009 


0.99 




Ss averaged in G' 


0.8768 


0.5584 


0.1107 


0.0412 


0.034 




Sg averaged in 5" 


1.119 


0.7287 


0.4493 


0.4751 


0.4448 




Ss averaged in S' 


1.062 


1.202 


1.009 


0.9293 


0.9002 


T2 


Sg averaged in G' 


1.031 


1.033 


1.01 


1.019 


0.997 




Ss averaged in G' 


0.823 


0.5146 


0.0842 


0.0388 


0.0377 




Sg averaged in 5" 


0.911 


0.2632 


0.0274 


0.0344 


0.0041 




Ss averaged in S' 


1.158 


1.201 


1.072 


1.049 


0.992 


T3 


Sg averaged in G' 


1.03 


1.057 


0.9868 


1.019 


1.009 




Ss averaged in G' 


0.869 


0.6141 


0.1164 


0.0085 


0.0067 




Sg averaged in S' 


1.001 


0.428 


0.0247 


0.0034 


0.0003 




Ss averaged in S' 


1.101 


1.177 


0.947 


0.8359 


0.783 



itself is consistent with unity. There seems a tendency of slightly lower values in the high 
luminosity especially for the mSUGRA case in T3. This may be due to varying errors of 
the significances by the statistical uncertainty of the Monte Carlo simulation. We checked 
this tendency vanishes if we exclude this uncertainty in the error of the significances. (See 
Table [Hi) In contrast, sm calculated in the wrong model decreases with higher integrated 
luminosity, indicating that more model points are discriminated from the other model. 

This fact is more clear in Table [TJ where the number of falsified model points using the 
criterion in Eq. f[T^ is presented. We see a large portion of the model points are falsified 
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Table 7: The number of falsified model points with respect to the integrated luminosity 
L. Eq. f lT^ is used for the falsification criteria. 





T 

Li 


i ID 


o ID 


1 n fK— 1 
iU ID 


oU ID 


1 nn fK— 1 
iUU ID 


Tl 


G pts (/1091) falsified m G 





1 


2 


1 


4 




G' pts (/1091) falsified in S 





323 


994 


1051 


1058 




S' pts (/754) falsified in G 





126 


199 


176 


169 




S' pts (/754) falsified in S 











5 


6 


T2 


G' pts (/1091) falsified in G 





1 


2 










G' pts (/1091) falsified in S 





394 


1049 


1089 


1089 




S' pts (/754) falsified in G 


10 


523 


753 


753 


751 




S' pts (/754) falsified in S 

















T3 


G' pts (/1091) falsified in G 








1 










G' pts (/1091) falsified in S 





136 


921 


1091 


1091 




S' pts (/754) falsified in G 





336 


752 


753 


754 




S' pts (/754) falsified in S 











1 


1 



from the other model at higher luminosity while most of them remain consistent within 
the home model. Let us see each tests in more detail. In Tl, many mSUGRA model points 
remain consistent with being a mGMSB model even at 100 fb"\ This is because without 
using the tau-jet mode, the slepton-NLSP GMSB model with nearly degenerated stau 
and the lightest neutralino may not be discriminated from mSUGRA models as described 
in the previous section. These mGMSB model points conversely remain unfalsified as 
a mSUGRA model. These models may be discriminated by using all five modes as in 
T3. The result shows all mGMSB model points are falsified from being in the mSUGRA 
model set at 30 fb"\ The T2 using combined significances also shows a good separation of 
models. In particular, it shows a better falsifiability in lower luminosity. This is because 
T2 uses only three variables for the test and has an advantage in statistics. In Table El 
we list the same result but without the Monte Carlo statistical error in the significances 
for T3. An important conclusion of this result is that even with the inclusion of the SM 
background and the effect of the detector efficiencies, the model discrimination among 
the GMSB and SUGRA models using the method gives a very good result at the 30 fb~^, 
i.e. at the time of the 5(j (or more) discovery. 
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Table 8: Same as Tables [6] and [7] for T3, without the Monte Carlo statistical error included 
in the significances. 





L 


1 fb-^ 


3 fb-^ 


10 fb"^ 


30 fb-^ 


100 fb^^ 


T3 


sg averaged in G' 


1.031 


1.053 


0.974 


1.009 


1.012 




ss averaged in G' 


0.856 


0.631 


0.093 


0.007 







Sg averaged in S' 


0.996 


0.394 


0.014 










Ss averaged in S' 


1.153 


1.277 


1.06 


1.041 


0.997 




G' pts (/1091) falsified in G 








1 










G' pts (/1091) falsified in S 





137 


959 


1091 


1091 




S' pts (/754) falsified in G 





360 


752 


754 


754 




S' pts (/754) falsified in S 


















In the above result, the model points are generated on a fiat grid in the parameter space 
and the tested models G' and 5" are subset of G and 5, respectively. There is a worry 
that the number of mis-falsified models is underestimated and there might exist small 
characteristic regions dropped out from the grid. To check this problem, we generate 
a new set of test model points randomly in the parameter space. The used regions in 
the parameter space is -100 < < 400 (TeV), 100 < Ag < 400 (TeV) and 5 < 
tan/3 < 50 for the mGMSB model and < mo < 1200 (GeV), 100 < mi/2 < 600 (GeV), 
5 < tan/3 < 50 and — 1 < Ao/ruQ < 1 for the mSUGRA model. After requiring the 
experimental constraints and the discovery condition {5a at 30 fb~^), we generated 885 
model points (denoted by G") for the mGMSB model and 1131 model points (denoted 
by S") for the mSUGRA model. We present the result of the discrimination test T3 in 
Table M The result is essentially same as the previous result for G' and S' and showing 



Table 9: The number of falsified model points by the significance test T3 for randomly 
chosen model points in G" and S". 





L 


1 fb-i 


3 fb-i 


10 fb-i 


30 fb-^ 


100 fb-i 


T3 


G" pts (/885) falsified in G 





3 


4 


4 


5 




G" pts (/885) falsified in S 





133 


809 


885 


885 




S" pts (/1131) falsified in G 





507 


1129 


1131 


1131 




S" pts (/1131) falsified in S 
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the underestimation is very small. 

5 Discussion and conclusion 

In this paper, we discussed a strategy of model discrimination between the mSUGRA 
and the mGMSB models at a very early stage of the LHC by comparing significances in 
multiple search modes. We defined the difference of the LHC signatures among model 
points in a quantitative form and show that the discrimination is possible with a relatively 
low integrated luminosity. Our estimation of the significances is conservative and the 
discrimination may be better in the actual experiments. This model discrimination is 
very important for both model-building and cosmology, since they predict completely 
different cosmological history and the dark matter. We used a relatively small number of 
modes in this paper. Although such a coarse graining of LHC signals gives only coarse 
information, it can discriminate the models qualitatively and guides to research at a much 
higher integrated luminosity. Though in this paper, we focused on only discrimination 
between the mSUGRA and mGMSB models, our method can be applied for more generic 
cases. For example, we can also discuss the discrimination within the mSUGRA models. 
The so-called co-annihlation and focus regions, which are favored from the viewpoint of 
dark matter, can be separated by using similar method. This method makes it possible 
to compare the LHC signals of generic BSM models qualitatively. 

Acknowledgement 

We would like to thank K. Hamaguchi and T. T. Yanagida for useful discussions. The 
work of E.N. is supported in part by JSPS Research Fellowships for Young Scientists. 

A Supplemental data of example model points 

In this appendix, we collect supplemental data of our example model points in Tab. Hlfor 
convenience. In Tab. [TO]and[TTl we list up the mass spectrum of the superparticles of the 
example model points together with their dominant decay channels. 

In Tab. [121 ^6 list the number of signal events for the example model points and the 
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Table 10: Superparticles' masses (in unit of GeV) and their dominant decay channels for 
the example model points PT 1 to 6 (mGMSB model) in Tab. |H From top to bottom, 
the particles' label represent gluino, down-type squarks, stop, sbottom, selectron, stau, 
charginos and neutralinos. The gravitino is denoted by G. 





PT 1 


Decay mode 


PT 2 


Decay mode 


PT 3 


Decay mode 


9 


1530.7 


xtm', X2W, Xim 


1530.1 


xfQQ', xIqq, Xim 


1319.4 


qq 


cLl 


2546.9 


gd.xiu 


2546.3 


gd, XiU 


1306.4 


XlU,X2U 


dR 


2414.7 


gd 


2414.0 


gd 


1252.4 


X\d 




2472.2 


gt, xtb 


2441.1 


gt, xtb 


1277.2 


Xlt,xlt,xtb 


h 


2208.4 


gt, xjb 


2210.9 


~ , —4-7 

gt, xjb 


1147.8 


Xl 0, X¥, xtb 


h 


2454.8 


gb, 


2428.0 


gb, x^t 


1261.0 


X2t,Xit 




2409.2 


gb 


2344.6 


gb 


1247.0 


&,& 




851.1 


— ~n ~n 

Xl X2e, Xie 


850.7 


Xl i", xie, Xie 


394.4 


~n 
Xie 




415.7 


Xie 


416.9 


Xie 


188.7 


Ge 




848.2 


Xl X2^, XV 


841.8 


Xl i", X2^, Xi^ 


392.1 


xy 




422.0 


Xir 


367.4 


Xir 


191.3 


Gt 


~+ 


877.7 


XihjZ, xW + 


860.8 


~ -4~ 1 / rz ~nTT7'-l- 

Xth/Z, xW + 


537.7 




xt 


529.8 


xW^ 


531.3 


TU, XiW+ 


423.8 


T T 7"-!- ~- /? 


xl 


877.3 


XjW ,X2^ 


860.2 


Xt^ ^Xlh 


539.8 


ii, vv 


xl 


867.6 


xtw-,xlz 


852.0 


xtw-,xlz 


483.0 


x'iz,ii 


xl 


529.5 


Xih 


530.8 


TT, Xih 


422.2 


Xih,ii 


X? 


273.4 


G-f/Z 


273.8 


Gi/Z 


240.6 


iRi,fir 




PT 4 


Decay mode 


PT 5 


Decay mode 


PT 6 


Decay mode 


9 


941.5 


xIqq, x\m, xtqq' 


1770.1 


qq 


1595.6 


qq 


dL 


1192.3 


gd,xlu 


1253.0 


X2U,xld 


1262.1 


X2U^ Xld 


dR 


1137.5 


gd 


1218.1 


xtid 


1220.2 


xtid 




1146.5 


xtb,X%xlt 


1242.6 


Xlt,xlt,xtb 


1230.4 


xtb,X%xlt 


ii 


1042.6 


Xib, X¥ 


1120.1 


xtb, Xlt, Xlt 


1122.6 


xtb, xlt, xlt 


h 


1139.4 


X2t^gb 


1220.5 


&,& 


1218.6 


&,& 


bi 


1093.5 


gb,xit 


1209.5 


&,& 


1178.0 


Xit,xlb 


ei 


378.7 


Xi^, X?e,X2e 


327.3 


fire 


350.1 


xie 


eR 


184.9 


X?e 


146.7 


Ge 


163.0 


TTe 


r2 


383.5 


riZ/Kx\r 


322.3 


fiZ/h 


354.9 


fiZ/h 


h 


137.5 


Gt 


149.5 


Gt 


114.1 


Gt 


xt 


473.0 


xlw+.xtz 


688.2 




616.8 


iu, pi 


xt 


304.8 


TV 


445.2 


i)i,xlw+ 


441.5 


fi U, VT 


xl 


472.9 


xtW-^xlh 


695.3 


z>z/, hi 


623.1 


ii, VV 


xl 


453.7 


xtw-,hT,xlz 


454.0 


X'lZ,hT 


454.6 


fiT 


xl 


304.1 


TiT 


447.3 




439.4 


fiT 


X? 


160.9 


TiT 


342.8 


fir,iRi 


306.3 


hr,iRi 
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Table 11: Same as Tab. [TOl but for example model points PT 7 to 10 (mSUGRA model) 
in Tab. |H Since the list of decay modes is very similar for PT 9 and PT 10, their lists 
are gathered. 





PT 7 


Decay mode 


PT 8 


Decay mode 


PT 9 


PT 10 


Decay mode 


9 


1157.2 


qq 


1156.7 


qq 


689.7 


689.1 


xtQQ', xlm, Xim 


dL 


1086.1 


XiU, xld 


1086.8 


XiU,X°2d 


1223.3 


1224.2 


gd, XiU 


d-R 


1043.4 


X\d 


1044.0 


xld 


1217.3 


1218.0 


gd 


i2 


1035.7 


xtb, 


1008.4 


xtb, xlt 


1031.5 


961.8 


gt, xtb 


ii 


813.5 


xtb, x\t, xtb 


816.1 


xtb, Xit, xtb 


755.4 


761.9 


Xtb,xlt 


h 


1038.5 


X2t, x\b 


1000.3 


X2t 


1204.5 


1076.4 


gb 


bi 


994.7 


Xit^&^xlb 


940.0 


Xlt,X2t,X^b 


1016.5 


941.1 


gb, & 




427.9 


Xie^Xii" 


427.9 


Xie^Xii" 


1111.7 


1112.0 


Xi^, xl(i 


en 


321.2 




321.6 


X?e 


1103.7 


1103.7 






426.9 


Xi-r, Xii^ 


423.9 


riZ/h,x'ir 


1107.5 


1027.9 




fl 


316.3 


X\r 


213.6 


Xtir 


1092.4 


920.4 


X\r 


xt 


650.5 


xlW^,xtZlh 


637.0 


X^W+,xtZ/h 


438.0 


404.0 


x^w+,xtz/h 


xt 


388.4 




390.6 




200.4 


201.1 


XiW+ 


xl 


651.0 


xtW'^xlh 


636.8 


xtW-,xlh 


438.1 


403.4 


xtW-,xlh 


xl 


636.4 


xtw-,xlz 


624.2 


xtw-,xlz 


422.3 


388.4 


xtw-,xlz 


xl 


387.9 


X\h, fiT 


390.0 


TiT 


199.8 


200.4 


xtiz 


X? 


205.7 


stable 


206.3 


stable 


105.8 


106.4 


stable 
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number of the background events for the search modes in Tab. [3l For the four jets with 
no lepton mode, the showed resuhs are event numbers under the hardest cut in Tab. [3l 
which optimizes the significance for most of the model points. 

Table 12: The number of signal events of the example model points for the five modes 
in Tab. |3] together with our estimation of the SM background. The event numbers are 
normalized for an integrated luminosity of 30 fb~^ at the 7 TeV run. 





0£ + 2 jets 


Oi + A jets 


27 


2r 


SS 2i 


PT 1 


0.3 


0.5 


45.3 








PT 2 





0.3 


44.9 


0.17 





PT 3 


9.1 


6.0 





6.04 


63.5 


PT 4 


67.4 


85.3 





65.1 


10.1 


PT 5 





10.7 





10.7 


57.8 


PT 6 


13.5 


9.0 





13.5 


57.5 


PT 7 


310.4 


79.0 


1.1 


1.14 


10.3 


PT 8 


333.8 


55.8 





2.37 


3.55 


PT 9 


71.1 


94.0 





2.30 


34.4 


PT 10 


51.84 


78.9 





9.01 


37.2 


SM BKG 


195.1 


44.9 


7.68 


14.8 


20.6 
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